To appear in 

Dynamics of Continuous, Discrete and Impulsive Systems 
http:monotone.uwaterloo.ca/ ^journal 



NONLINEAR DYNAMICS IN QUANTUM 
PHYSICS-THE STUDY OF QUANTUM CHAOS 

L.A. Caron a , H. Kr6ger a , X.Q. Luo 6 ' c , 
G. Melkonyan ' 1 and K.J.M. Moriarty d 

a Departement de Physique, Universite Laval, Quebec, Quebec G1K 7P4, Canada 
b CCAST (World Laboratory), P.O. Box 8730, Bejing 100080, China 
c Departement of Physics, Zhongshan University, Guangzhou 510275, China 
^Department of Mathematics, Statistics and Computer Science, 

Dalhousie University, Halifax, N.S. B3H 3J5, Canada 
1 Talk given at DCDIS conference in London (Canada) by G. G. Melkonyan, Email: 

gmelkony @phy. ulaval . ca 

Abstract: 

We use the quantum action to study the dynamics of quantum system at fi- 
nite temperature. We construct the quantum action non-perturbatively and 
find temperature dependent action parameters. Here we apply the quan- 
tum action to study quantum chaos. We present a numerical study of a 
classically chaotic 2-D Hamiltonian system - harmonic oscillators with an- 
harmonic coupling. We compare Poincare sections for the quantum action 
at finite temperature with those of classical action. 
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1 Introduction 

Chaotic phenomena are found in microscopic physics ruled by quantum me- 
chanics - an area presently under very active investigation. For reviews see 
Refs.jjJ Jl2|, (l3[ 0. Examples are the hydrogen atom in a strong mag- 
netic field [p~o[| and the quantum mechanical stadium billiard j23|. Billiard 
like boundary conditions have been realized experimentally in mesoscopic 
quantum systems, like quantum dots and quantum corrals, formed by atoms 
in semi-conductors |u|. 

Classical chaos theory can not simply be taken over to quantum physics 
(due to Heisenberg's uncertainty relation). Hence workers have tried to char- 
acterize quantum chaos in ways alternative to classical chaos. The following 
approaches to quantum chaos have been considered: 

(i) Gutzwiller's trace formula. Gutzwiller has established a relation 
between the density of states of the quantum system and a sum over classical 
periodic orbits (periodic orbit quantisation). The trace formula has been 
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applied in the semi-classical regime (e.g. highly excited states of atom). 
Wintgen Q applied it to the diamagnetic hydrogen system and extracted 
periodic orbit information from experimental level densities. 

(ii) Effective action. Another possibility is to use the conventional effec- 
tive action [|l9[ |(|. An effective action exists also at finite temperature 
The effective action has been introduced in such a way that it gives an ex- 
pectation value <(/)>= (j) class which corresponds to the classical trajectory 
and which minimizes the potential energy (effective potential). Thus one 
can obtain the ground state energy of the quantum system from its effective 
potential. Because the effective action has a mathematical structure similar 
to the classical action, it looks like the ideal way to bridge the gap from 
quantum to classical physics and eventually solve the quantum chaos and 
quantum instanton problem. Cametti et al.jij have computed in Q.M. the 
effective action perturbatively via loop (ft) expansion and studied quantum 
chaos in the 2-D anharmonic oscillator. Perturbation theory forces to termi- 
nate such series at a finite (low) order. The difficulty lies in the estimation 
of the remainder. This poses a problem in the context of chaos, because 
chaotic systems are known to exhibit sensitivity not only to initial conditions 
but also to parameters of the system. There are higher loop corrections to 
the effective potential V e ^ as well as to the mass renormalization Z. The 
kinetic term of the effective action is given by an asymptotic series of higher 
time derivatives (infinite series of increasing order). The problem in inter- 
preting T as effective action is that those higher time derivatives require more 
initial/boundary conditions than the classical action. 

(hi) Bohm's interpretation of quantum mechanics. An interesting ap- 
proach has been taken by Iacomelli and Pettini . They have used Bohm's 
H interpretation of quantum mechanics, which expresses the Schrodinger 
equation in polar form. The radial part satisfies a continuity equation of a 
velocity field and the phase satisfies a Hamilton- Jacobi equation. The "tra- 
jectories" of the quantum system are equivalent to Lagrangian trajectories 
of a fluid velocity field. This allows to apply the concept of classical chaos 
theory and in particular to compute Lyapunov exponents. In Ref.[|l5| this 
framework has been used to study the chaotic behavior of the hydrogen atom 
in an external magnetic field and different regimes of quantum chaoticity have 
been observed. 

(iv) A conjecture by Bohigas et al. J|] states that the signature of a 
classical chaotic system, described by random matrix theory, is a spectral 
density following a Wigner distribution. 

Strictly speaking, there is no universally accepted definition of quantum 
chaos. Many roads have been travelled searching for "irregular" signals in 
quantum systems. In the quest for quantum chaos, it has been popular 
to examine systems which exhibit classical chaos. However, Jona-Lasinio 
et al. [^0| found chaotic behavior in a quantum system without classical 
counterpart, i.e., a quantum many-body system undergoing multiple resonant 
tunneling. Jona-Lasinio and Presilla p2| further studied quantum many- 
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Figure 1: Path integral takes into account all possible paths between two 
boundary points xi n , Xfi 



body systems of interacting bosons in the thermodynamic limit. They found 
chaotic signals when the interaction strength drives the system away from 
integrability. 

However, the study of quantum chaos lacks from the following shortcom- 
ings: (i) In quantum chaos, we have no "local" information of the degree 
of chaoticity, being available in classical chaos via Lyapunov exponents and 
Poincare sections from phase space, (ii) Also little is known about the role 
of temperature in quantum chaos. For example, the analysis of level densities 
is insensitive to temperature. To overcome those problems, we adopt the 
strategy of building a bridge between quantum physics and classical physics. 
This means a relation involving the quantum transition amplitude and the 
classical action. 

In the section || we present the definition of the quantum action ]l6| . In 
the section || we describe methods of construction of the quantum action. 
The dynamics in 2-D anharmonic oscilator by use of the quantum action is 
considered in the section 0. 



2 Quantum action 

Feynman's path integral formulation of quantum mechanics j|, |25| is helpful 
in introducing the ideas of the quantum action. In this formulation the 
quantum mechanical transition amplitude or Greens function is defined as 

G(xfi,tfi,Xin,Un) = J d[x] exp^S[x]j, (1) 

where S is classical action, a^ n , Xfi are initial and final boundary points, 
Umtfi are initial and final times. The integral is understood as an integral 
over all possibles paths connecting initial and final boundary points (see 
Figure [l|). In some exceptional cases (harmonic oscillator and free particle) 
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integral (Q) can be presented as 

G(xfi,tfi,Xi n ,t in ) = Z exp (j^ s l x ci]j 



Xfi,t fi 

(2) 

Xin ,tin 



where x c \ is the classical path which minimizes the classical action and solves 
Euler-Lagrange equations of motion. The definition of the quantum action 

is guided by Equation (|2|). 
Conjecture: For a given classical action S = J dt^x 2 — V(x) there is a 
quantum action S = J dt^x 2 — V(x), which allows to express the quantum 
mechanical transition amplitude G as 



G{xfi,tfi\Xi n ,t in ) = Zexp ( ^S[x c i] 



X fi ,tfi 



(3) 



Here x c \ denotes the classical path, such that the action S[x c i) is minimal 
(we exclude the occurrence of conjugate points or caustics). Z denotes the 
normalization factor corresponding to S. Eq.(^) is valid with the same action 
S for all sets of boundary points Xfi, Xi n for a given time interval T = 
tfi — Un- The parameters of the quantum action depend on the time T. Any 
dependence on Xfi,xi n enters via the trajectory x c \. Z depends on the action 
parameters and T, but not on Xfi,Xi n . 

This action has a mathematical structure like the classical action, and 
takes into account quantum effects (quantum fluctuations) via renormalized 
action parameters. The quantum action has the following attractive prop- 
erties [ p?r| : (i) It can be computed non-perturbatively, via numerical simu- 
lations, (ii) It is defined for arbitrary finite temperature, (iii) In the zero- 
temperature limit analytic expressions exist for the quantum action, provided 
that the classical potential is sufficiently symmetric (parity symmetric in D=l 
and rotationally symmetric in D=2,3). (iv) In a numerical study of the 1- 
D anharmonic oscillator, we searched for the presence of higher order time 
derivative terms in the kinetic term of the quantum action. No signal for 
such term was found [|l7| , in contrast to the effective action formalism, where 
higher derivative terms persist. The quantum action has proven useful for a 
rigorous definition and quantitative computation of quantum instantons . 

The quantum action at finite time T can be interpreted as action at finite 
temperature. The partition function in quantum mechanics requires to go 
over from real time to imaginary time (t = —it) and to impose periodic 
boundary conditions In Ref.[[l7j we have shown that the expectation 

value of a quantum mechanical observable O at thermodymical equilibrium 
can be expressed in terms of the quantum action along its classical trajectory 
from x in , Pin = to x, (3, given by 
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Xcl(Pin) = Xin, X c i(P) = X , (4) 

where T denotes the absolute value of imaginary time. The parameters /?, 
temperature r and T are related by (3 = ^7 = T/7i. Numerical studies with 
the quartic potential in 1-D have shown that the parameters of the quantum 
action are temperature dependent [|l6[ |l8], |l7[ . 



3 Methods of construction of the quantum ac- 
tion 

3.1 Analytic form of quantum action at zero tempera- 
ture 

From the experimental point of view, the ground state properties of a quan- 
tum system are often much better known than excited states. Zero tempera- 
ture corresponds to large time and from the point of view of chaos theory, it 
is interesting to study dynamical behavior at large times. Examples: (a) Par- 
ticles moving in storage rings, (b) In classical Hamiltonian chaos, there are 
quasi-periodic trajectories, remaining for a long time in the vicinity of regu- 
lar islands. What is the corresponding behavior in quantum mechanics? (c) 
What is the quantum analogue of the KAM theorem? (d) The construction 
of Poincare sections and Lyapunov exponents involves long time trajectories. 
Likewise the quantum analogue of Poincare sections and Lyapunov exponents 
will require long time evolutions. 

In the following we work in imaginary time. In order to avoid notational 
confusion, we continue to use the time parameter T. In particular, T — > 
denotes the high temperature limit, where the quantum action approaches its 
classical counter part. T — > 00 denotes the zero-temperature limit, where the 
quantum action is determined by the Feynman-Kac formula. In Ref.JlTj the 
following analytic transformation law at temperature zero has been derived 
in 1-D relating the classical action S = J dt\mx 2 + V(x) to the quantum 
action S = J dt\mx 2 +V(x), 

2m(V(x) - E gr ) = 2m(V(x) - v~o) - ^ 2m(F(a?) ~ ^ Lgn(x). (5) 

y/2m(Y(x) - v~o) 

This differential equation relates the classical action to the quantum action at 
zero temperature limit. Moreover, the following relation between the ground 
state wave function and the quantum action holds 

^gr(x) = — e Jo v E gr =v , (6) 
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Eqs.(||||) are valid under the assumption of a parity symmetric potential 
V{x) with a unique minimum at x = 0, V"(0) = %. Similar laws hold for 
rotationally symmetric potentials in D=2 and 3. For a number of systems, 
where the ground state energy and wave function is known analytically, it 
has been shown |l7j that the quantum action reproduces exactly the ground 
state energy and wave function. Examples in 1-D are: Monomial potentials 
and the inverse square potential. Examples for radial motion in 3-D are: 
Monomial potentials and Coulomb potential (hydrogen atom) in the sector 
of non-zero angular momentum. In particular, the ground state energy co- 
incides with the minimum of the quantum potential V and the location of 
the minimum of the quantum potential coincides with the location of the 
maximum of the wave function. One should note that the transformation 
law, Eq.(||), determines mass x potential, but not both, mass and potential 
individually. The determination of the renormalized mass requires a non- 
perturbative numerical calculation. 

3.2 Calculation of quantum action by fit to the quantum 
mechanical transition amplitude 

A method to obtain the quantum action follows directly from the definition 
of the quantum action. Equation (||) must hold for all boundary points for 
a given time interval T = tfi — U n . A parametrization of the quantum 
action as a polynomial by taking into account symmetry of the system is 
very useful. It permits to fit the quantum action to the quantum mechanical 
transition amplitude. After parametrization one minimizes the sum of square 
differences 



to calculate the parameters of the quantum action. G(xfi,T, o^ n ,0) is the 
exact quantum mechanical transition amplitude obtained from the solution 
of the Schrodinger equation. 

3.3 Calculation of quantum action by renormalization 
group flux equation 

In quantum field theory action parameters depend on cut-off parameter A or 
on spatial and time lattice spacing. This dependence is governed by a renor- 
malization group equation. In constructing the quantum action we are close 
to the continuum limit. As scale and temperature dependence of the action 
are similar, one can write an equation which governs the action parameters 
when temperature is varied [0. Following our conjecture this equation 
reads 
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(7) 



X in ,Xfi 




(8) 
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E is given by the quantum action along its classical trajectory and (3 = T/h 
as before. 



~ ~ f 171 • 2 ~ 

E/3 = Sp[xd] = / dt—^x cl + V(x c i) 
Jo Zri 



(9) 



By parametrizing the quantum potential in terms of Vk\x\ k correspond- 
ing to the symmetry of V(x) from equation (||) we have a system of linear 
differential equations for the quantum action parameters Vk,fh and Z. 

= E/3 [m(J3), %(J3), V! (13), . . . , x, [3] . (10) 

One can solve this system of differential equations either by giving initial 
conditions at (3 = (classical limit) or at (3 = oo. Like equation (||) this sys- 
tem of differential equations connects the classical and quantum potentials. 
It determines the flow of the action parameters as temperature is lowered or 
increased. 



3.4 Numerical results 

First we will try to estimate the precision of numerical calculations of the 
quantum action in the zero-temperature limit (long time limit) by contrasting 
it with the analytic expressions of Eqs.(||J^). We consider in 1-D the classical 
action 

S = J dtfx 2 + V(x), V(x) = v 2 x 2 + v 4 x 4 , 

m = 1, v 2 = 1, va = 0.01 . (11) 

The Q.M. Hamiltonian has the following ground state energy, 

E gr = 0.710811 . (12) 

For the quantum action we make an ansatz 

S = J dt^x 2 + V(x), V(x) = vo + v 2 x 2 + V4X 4 + v 6 x 6 . (13) 

From a fit of the quantum action to Q.M. transition amplitudes (expres- 
sion (0), for more details see Ref.Jl^]) we find the action parameters at 
T = 4.5 

rh = 0.9990(7), % = 0.79863(1), v 2 = 1.013(2), 
v 4 = 0.0099(9), v 6 = 0.0000(3) . 

The parameters m, v 2l £4, vq have asymptotically stabilized at about T = 4.5. 
However, the parameter £0 has a slow fall-off behavior for T — » 00. For large 
T we have fitted £0 to the numerical data by 

v (T)^ T ^ oo A+^ + ^ , (14) 
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Figure 2: 1-D anharmonic oscillator. Differences between ground state wave 
functions obtained from Schrodinger equation and from quantum action. 



and find A = 0.710819. Thus the asymptotic extrapolation vo(T — > oo) = 
0.710819 is very close to E gr , given by Eq.(|T^). This confirms E gr = vq, pre- 
dicted by Eq.(^|), i.e. the ground state energy is given by the minimum of the 
quantum potential. Let us further check the consistency of other parameters 
of the quantum potential. This can be done by using the transformation 
law Eq.(|5|). Because this transformation law does not determine m, we have 
taken rh from the numerial data. The transformation law then determines 
the parameters of the quantum potential. We find, 

v 2 = —\mE Qr /h] 2 = 1.01151 
m 



v 4 = l^^[v 2 -v 2 ^] = 0.009967 
3 n m 



v 6 = 1~ T — T 1^4 



t . ■ r t ^. ^4^]=2.89 10- 7 . (15) 
4 V2 5 a m 

This is consistent with the numerical results, given by Eq.(|lJ). Finally, let 
us look at the ground state wave function. As reference solution, we have 
solved the Schrodinger equation with the original Hamiltonian. Secondly, we 
have computed the wave function from the quantum action via Eq.@, using 
the numerical data given by Eq. (|l4|) . The comparison is given in Figure ||, 
showing good agreement. 

Now let us consider the quantum action at finite temperature. The re- 
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Figure 3: Differences between quantum action parameters and effective ac- 
tion. Potential V = 0.5x 2 + Ax 4 . 



suits obtained by minimizing the error in the expression (|7|) (direct fit) show 
that at high temperature parameters of quantum and classical actions co- 
incide and at low temperature the parameters of the quantum action are 
saturated ]l6| . Figure || presents a comparison of the quantum action with 
the conventional effective action. We see at small coupling that the quantum 
and effective action are similar but at high values of coupling there is some 
discrepancy. Finally, we presnt a comparison of the methods to determine 
the quantum action- via fit (Q) and via solution of the renormalisation group 
flux equation (JsJ) . Such comparison is presented in Figure |[ Two different 
methods reproduce the same result. 

In summary of this section, we have shown for the zero-temperature limit 
(regime of the Feynman-Kac formula) that the transition amplitude can be 
parametrized by the quantum action. A comparison of analytical versus 
numerical computation of the quantum action has shown good agreement. 
We have shown that in effective and quantum actions can coincide at small 
coupling constants. Also we show that solutions obtained by two different 
methods (direct fit to quantum mechanical transition amplitude and flux 
equation) agree very well. 

4 Quantum chaos in 2-D anharmonic oscillator 

As is well known 1-D conservative systems with a time-independent Hamilto- 
nian are integrable and do not display classical chaos. An interesting candi- 
date system to display chaos is the K-system, corresponding to the potential 
V = x 2 y 2 . It describes a 2-D Hamiltonian system, being almost globally 
chaotic, having small islands of stability 0. However, according to a sugges- 



9 



I ,004 

1 ,25 
1 ,00 

0,7^5 - 



:!!!!! 



1 1 1 



0,50 



0,5 1 ,0 1,5 2,0 



,5 3,0 3,5 4 



V 4 

0,025 
0,020 
0,01 5 
0,01 O 

m 
0,998 
0,996 
0,994 
0,992 

0,990 
O 



Figure 4: Quantum action parameters obtained from the direct fit and from 
solution of the flux equation. Potential V = 0.5x 2 + O.Ola; 4 . Lines: flux 
equation and symbols: fit. 



tion by Pullen and Edmonds it is more convenient from the numerical 
point of view to consider the following system, which also exhibits classical 
chaos. It is a 2-D anharmonic oscillator, defined by the following classical 
action, 

S = Io dt ^rn(x 2 +y 2 ) + V(x,y), V(x, y) = v 2 (x 2 + y 2 ) + v 22 x 2 y 2 

rn = 1, v 2 = 0.5, v 22 = 0.05 . (16) 

We work in imaginary time and the convention Ti = ks = 1 is used. For 
the corresponding quantum action, we make the following ansatz, compat- 
ible with time-reversal symmetry, parity conservation and symmetry under 
exchange x <-> 

f T 1 

S = dt -m(x 2 +y 2 )+V(x,y), 

V = vo + v 2 (x 2 +y 2 ) +v 22 x 2 y 2 + v 4 (x 4 + y 4 ) . (17) 

We have determined numerically the parameters of the quantum action for 
transition times from T = up to T = 4, where a regime of asymptotic 
stability is reached. For T = 4 corresponding to temperature r = 0.25, we 
found 

m = 0.99814(6), v = 1.29373(2), v 2 = 0.5154(3) 

v 22 = 0.0497(1), v 4 = -.0008(7) . (18) 
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We have also included terms xy, xy, xy 3 + x 3 y. Their coefficients were found 
to be very small (compared to machine precision). Also we included terms 
x 2 y 4 + x A y 2 , x 4 y 4 . Those coemcients were found compatible with zero within 
error bars. The quantum action slightly modifies the parameters V2 and the 
parameter £22- 

Now let us look at the chaotic behavior of the classical and the quantum 
action. We have computed Poincare sections from the classical action and 
from the quantum action for a variety of temperatures. The equations of 
motion of the quantum action have been solved using a 4-th order Runge- 
Kutta algorithm, and Henon's algoritm was used to compute the Poincare 
sections. Let us discuss what we expect. First, the classical system under 
consideration is chaotic due to the anharmonic interaction term. This term 
is small compared to the harmonic oscillator term. This property proliferates 
to the quantum potential. Here "smallness" is meant in the following sense: 
The quantum system at zero temperature is determined by the ground state 
properties. We have computed ro = (r) gr , the radius of the ground state 
wave function (corresponding to the Bohr radius ao in the hydrogen atom). 
We obtain ro = 0.8766. It turns out that the anharmonic term is smaller 
than the harmonic term for any length r < tq. Thus one expects roughly 
the same "amount" of chaotic behavior in the quantum system as in the 
classical system. Second, renormalization (i.e. change of action parame- 
ters when going from the classical action to the quantum action) tends to 
increase harmonic terms, i.e. drive the system to a Gaussian fixed point. 
At a Gaussian fixed point the system is integrable. i.e. non-chaotic. This 
effect has been observed in a number of numerical studies in 1-D systems 

[l6|, [l8|. In a system with a symmetric double well potential this effect 
was found to lead to a softening of the double well shape of the quantum 
potential (lower barrier, more narrow minima, even vanishing of the double 
well) . This translates to a softening or even an " evaporation" of quantum in- 
stantons compared to the classical instantons [|l8| . This effect is present also 
in the system under consideration here. We find Sv2 = ^2 — ^2 = 0.0154(3) 
compared to ££22 = ^22 — ^22 = —.0003(1). We expect that this effect may 
lead to a slight softening of the chaotic behavior in the quantum system. 
Third, a weak anharmonic term means that the quantum potential reaches a 
regime of asymptotic stability as function of T already for small T (T « 4). 
That is, when going from high to low temperature, already at relatively high 
temperature the asymptotic regime (Feynman-Kac limit) is reached. 

We computed the Poincare sections for the quantum action at temper- 
ature r = 0.25 (corresponding to T = 4) for energies E = 10, 20, 50. One 
should keep in mind that the classical action at T = is equivalent to a 
quantum action at temperature r = 00. Poincare sections compared to their 
counter part from the classical action for energies E = 10, 50 are shown in 
Figures || and ^. One observes that the quantum system also displays chaos, 
and the Poincare sections are slightly different from those of the classical 
action. Like in the classical case also in the quantum system the "amount" 
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Figure 5: 2-D anharmonic oscillator. Poincare sections at energy E = 
Left part classical action, right part quantum action at temperature r = 




Figure 6: Same as Figure ||, but E = 50 
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of chaoticity increases with increasing energy. Moreover, the difference be- 
tween classical and quantum Poincare sections becomes more accentuated 
with increase of energy. 

5 Summary 

We have discussed the use of the quantum action, which can be considered 
as a renormalized classical action at finite temperature and methods of its 
construction. Based on the quantum action we suggest a proper definition of 
quantum chaos and explore it numerically. As example, we have considered 
in 2-D the harmonic oscillator with a weak anharmonic coupling and com- 
puted the quantum action at various temperatures. We compared Poincare 
sections of the quantum action with those of the classical action and found 
some difference in their chaotic behavior. 
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